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We use an adiabatic approximation in terms of instantaneous resonances to study the steady- 
state and time-dependent transport properties of interacting electrons in biased resonant tunneling 
heterostructures. This approach leads, in a natural way, to a transport model of large applicability 
, consisting of reservoirs coupled to regions where the system is described by a nonlinear Schrodinger 

■ equation. From the mathematical point of view, this work is non-rigorous but may offer some 

fresh and interesting problems involving semiclassical approximation, adiabatic theory, non-linear 
Schrodinger equations and dynamical systems. 
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I. INTRODUCTION 



Man-tailored semiconductor heterostructures Q offer, for the first time, the possibility to test quantum mechanics 
at a mesoscopic level ||^ . The scenario of systems which can be investigated is so rich that the art of their realization 
, deserves the name of quantum design. 
' In the simplest case, a quantum designer can grow sandwiches of different semiconductor alloys by choosing the 
\^ ' number of atomic layers for each kind of alloy. In the resulting heterostructure, the conduction band profile along the 
On growth direction forms steps whose height can be continuously varied by a proper choice of the alloy composition. 

' Typical widths and heights are of the order of tens of A and tenths of eV, respectively. 
^ . At low-temperature, the mean free path of carriers for scattering from crystal impurities is of the order of 10* A 
C ' and for heterostructures smaller than this size the electric transport along the growth direction is a phase coherent 
quantum scattering from the conduction band discontinuities [^|. Due to the translational invariance in the plane 
^ • orthogonal to the growth direction, the problem is one-dimensional. Moreover, the carriers are described by an 
Q ' effective mass which accounts for the microscopic scattering with the periodic crystal sites and their wave function is 
O • an envelope wave function Q|. 

^ ' In a homogeneous neutral conductor, the electron-electron interaction can be taken into account by a renormal- 
ization of the carrier effective masses |5|] and one deals with a transport problem like in a noninteracting case. In 
a heterostructure, even as simple as that described above, the breaking of translational invariance in the transport 
direction allows the electric neutrality to be locally violated. The corresponding interaction potential, obtained, at 
Hartree level, by solving a proper Poisson equation, can strongly modify the transport properties. The example of a 
double barrier heterostructure with the exterior regions doped with donors is illuminating [|6|. Due to tunneling, elec- 
trons populate the resonance(s) created by the double barrier and the region between the barriers becomes negatively 
charged. This generates an electric potential which decreases the tunneling probability of electrons in the double 
barrier region. As a consequence, current oscillations on the picosecond scale 1^,^ and chaotic behavior without 
classical counterpart |^ have been predicted in a ballistic configuration in which electrons are injected at some chosen 
energy. 

Experiments with ballistic electrons are difficult and measurements became available only recently | [To[ |. Techno- 
logically simpler is the case of biased heterostructures where transport is due to the presence of reservoirs at thermal 
equilibrium with different chemical potentials. Manifestations of the electron-electron interaction are known also in 
this configuration. For example, hysteresis in the current-voltage characteristics of double barrier heterostructures 



have been observed |11| and recognized as a consequence of the accumulation of electrons in the resonance [|ll|-[l5[. 
In this case, however, one has the theoretical problem of attaching reservoirs at thermal equilibrium to a piece of 
conductor where quantum coherent transport takes place. 

In the recent paper [T^ ] we proposed an approach to this problem based on a mathematical method earlier applied 
in the framework of ballistic transport |jl^ . We showed that for heterostructures with a single resonance our approach 
allows one i) to obtain steady-state voltage-current characteristics having hysteresis or not in agreement with the 
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experimental results (iSj and ii) to predict time-dependent properties analogous to those studied in optically bistable 
systems . Here, we develop the general mathematical scheme of this approach and discuss the case with several 
resonances where multistability phenomena can take place as in superlattices 1 20 2^ . 

For simplicity, consider the one-dimensional double barrier heterostructure discussed above. The idea is that due to 
the presence of resonances the corresponding Schrodinger problem can be divided in two parts: a Schrodinger equation 
for the barrier region and one for the exterior space, the two being weakly coupled by tunneling. This decomposition 
corresponds to the schematization of the transport process as a coherent process fed by reservoirs. In the exterior 
space (reservoirs) homogeneous and neutral, the electron-electron interaction is neglected and thermal equilibrium is 
taken into account by considering a continuous set of energy eigenstates distributed according to the Fermi statistics. 
In the barrier region (coherent conductor), the Coulomb interaction is included in a self-consistent potential obtained 
by solving the Poisson equation associated to the local charge density. Under the assumption that the barriers are 
wide enough, the corresponding nonlinear Schrodinger problem is discussed in two steps. In the first step we eliminate 
the potential well between the two barriers, by artificially increasing the potential there, and we solve the Schrodinger 
equation asymptotically for the new potential by means of WKB-expansions. The resulting solution is then very small 
near the (filled) potential well, so we get only a small error in the Schrodinger equation when we go back to the true 
potential. In the second step we correct for this small error by adding a wave function concentrated near the potential 
well. Assuming a priori that the charge in the well changes slowly with time, the correcting wave function can be 
expected to be large only at energies close to the resonances, and be well approximated by some linear combination 
of the resonant states. 

In most of the paper we discuss the case in which only one resonance participates. The validity of this one-mode 
approximation has been tested numerically with excellent results in the ballistic configuation of Here, the 

coefficient of the one-mode approximation obeys an ordinary differential equation with respect to time in the infinite 
dimensional space of square integrable functions of energy. We study the stationary points of the corresponding 
vector field and their nature, whether they are attractive or not, and arrive at quite neat answers. For solutions of the 
dynamical system which have existed as bounded solutions for a long time and in a suitable asymptotic limit (of wide 
barriers) we derive a simplified scalar differential equation for the evolution of the sheet density of electrons trapped in 
well which gives good global understanding of the more complete dynamical system. Using these results, we are able 
to discuss the phenomenon of hysteresis and we support and illustrate the discussion with several numerical results. 
The discussion includes the evolution of solutions away from fixed points which necessarily appears when there is 
hysteresis. We also discuss the case of several resonances, and get analogous results. 

From the mathematical point of view, the present paper could be a starting point for rigorous work on some fresh 
problems, involving semiclassical analysis, adiabatic theory, non-linear Schrodinger equations and dynamical systems. 
A strong motivation for such an enterprise is the fact that the theory of electric transport in semiconductor devices 
offers many problems similar to that one we illustrate here ||2^ . 

The plan of the paper is as follows. In Section |l^ we define the model. In Section IH we review the WKB expansion 
for slowly varying potentials. In Sections IV and ^ we determine the driving term and the ground resonant state, 
respectively, within the WKB approximation. The central equation of our paper is derive d in Section ^ and the 
general properties of the associated fixed points and linearizations are discussed in Section VII . In Section |VIII we 
introduce an approximation valid in the limit of small resonance width and discuss the corresponding fixed point 
solutions and linearizations. In section IX we obtain a simplified differential equation describing the dynamics of the 
electron density in the well. A qualitative discussion of the hysteresis phenomenon in comparison with numerical 
results is given in Section In Section XI wc finally consider the case with several resonances. 



II. DEFINITION OF THE MODEL 



Let us consider a heterostructure whose conduction band profile consists of two barriers of height Vq located in 
[a, b] and [c, d] 



(2.1) 



with a < b < c < d along the growth direction x. We wish to evaluate the transport properties of this device when 
a bias energy AV is applied between the emitter (x < a) and collector (x > d) regions uniformly doped. Due to 
doping, the band of conduction electrons formed in the emitter and collector regions is characterized by a Fermi 
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where n^i is the net donor concentration. We will use everywhere effective atomic units 
: 2 a~g^ , where m* is the electron effective mass and e the dielectric constant. In these 



units, every physical quantity is expressed in terms of the effective Bohr radius as = ^ e/{m*e ). Assuming an ideal 
heterostructure homogeneous in the plane yz parallel to the junctions (and orthogonal to the growth direction x), the 
single-electron momenta ky and are conserved quantities. As a consequence, the single-electron wavefunction at 
energy E + E^^, where E'y ~ ky + kl, can be factorized as (j){x, t, E) x(y, z, t, with 



x{y,z,t, Ell) 
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(2.2) 



We will assume periodic boundary conditions in a two-dimensional region A so that the momenta ky and k^ are 
quantized as in a real device having finite lateral area of size A. The time dependent Schrodinger equation for the 
single-electron wavefunction at energy E along the x direction is 



-idt -dl + V,i,{x) + U{^, x)] 4>{x, t, E) = 0, 



(2.3) 



where U{(j),x) takes into account the applied bias and, at Hartree level, the electron-electron interaction. Assuming 
ideal metallic behavior in the emitter and collector regions, i.e., neglecting the formation of accumulation and depiction 
layers, U (0, x) can be obtained as solution of the Poisson equation 



a2;7(0,x) = -8™sV(0) 



(2.4) 



with Dirichlet boundary conditions U {(f>, a) — and U {4>, d) — — AV^. The density p takes into account all the electrons 
in the occupied energy states and depends only on the wavefunction component 0. Indeed, if the emitter and collector 
regions are at thermal equilibrium with temperature T we have 



p = 2 / dEj2 \cbix,t,E) xiy,z,t,Eii)ni 



dEg{E) |0(x,t,£;)|2 



(2.5) 



where the factor 2 takes into account the spin degeneracy. Energies are me asured from the bottom of the emitter 
conduction band and the lower integration bound = in the first line of ( ^.5[ ) stems from the fact that for Ep <^ AF, 
as we will assume, only electrons from the e mit ter conduction band can penetrate the region [a, d] where the electron 
density is of interest. In the second line of (2^) this lower bound is absorbed in the definition of g{E) by a Heaviside 
function 9{E). The function g{E) can be explicitly evaluated by approximating the sum over the parallel degrees of 
freedom with an integral 



g{E) = e{E) 2 
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(2.6) 



Note that the chemical potential at temperature T in the Fermi function has been approximated with its value at 
T = 0, i.e., the Fermi energy determined by the net donor concentration. 

In general, the solution of (2.4) can not be handled analytically. We will suppose that, due to the accumulation of 
electrons in the well with sheet density 



s(</.) 



dE g{E) 



{c+d)/2 



(a+b)/2 



dx \(t>(x,t,E)f 



(2.7) 



ideal metallic behavior in the well \b, c] and ideal insulating behavior in the barriers [a, b] and [c, d] hold. This is 
equivalent to approximate (2.4) with 



dlU{(l}, x) = -87ra5\s(0) [B5{x - b) + Cd{x - c)] 



B + C=l 



(2., 



and the condition that dxU{(j), x) = for b < x < c. In this case U (0, x) becomes a piece- wise linear function of x with 
dxU{(j), x) having jump discontinuouities at x — b and x = c. The total potential Vch + U in ( |2.3| ) is better rewritten 
asV + W where 
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X < a 

Va- AV{x - a)/e a<x<b 

V{x) = { -AV{b- a)/e b<x<c (2.9) 

Vn- AV{b~a + x-c)/e c<x<d 
-AV x>d 

gives the band profile modified by the external bias and 

X < a 

{x — a){d — c)/£ a < X < b 

Wis,x) = 8TTa-^ s{(j)) { {b-a){d-c)/t b<x<c (2.10) 

{b-a)id-x)/£ c<x<d 

X > d 

depends on the wavefunction (p through the sheet density of electrons in the well s ((/)). Here £ — b — a + d — c. The 
potentials V{x) and W{s,x) are shown in Fig. 1. 

We will try to solve the nonlinear partial differential equation 

[~idt ~dl + V{x) + W{s, x)] (f>{x, t, E) = 0, (2.11) 



where s((^) is given by (2.7), in two steps. Let Vf^\\{x) = V{x) + Vol[6,c](2;) be the potential obtained by filling the well 



[6, c]. Here l\jr,^c]{x) is the characteristic function of the interval [6, c]. First we solve 

[~idt -dl + Vmi{x) + W{s, x)] fiix, t,E)^0 (2.12) 

and then we look for (f> in the form (j) ^ jl + i> where i> should solve 

[-idt -dl + V{x) + W{s,x)] i){x,t,E) = Vo l[b^c]{x)fl{x,t,E). (2.13) 

The wave function /2 describes an electron at energy E which is delocalized in the emitter and collector regions and 
has an exponentially small probability to be found in the forbidden region [a,d\. The wave function i) describ es the 



localization, driven by p,, of the same electron in the well [fe, c]. The wave function of the original problem (2.11) 
can be approximated by or /i inside or outside the two barriers, respectively, with an error which is exponentially 
small in the limit of wide barriers [ p^f 



To evaluate fl we will use a WKB approximation in the forbidden region [a,d]. Equation ( 2.13| ) will be treated 



with a one mode approximation in which D is assumed proportional to a resonant state corresponding to the potential 
V + W. To evaluate this resonant state and the corresponding resonance, we will again use a WKB approximation. 
In both cases, the justification of using a WKB approximation stems from the fact that Van + W and V + W are 
slowly varying potentials in the barriers regions if 6 — a and d — c are large while AV and s remain bounded. 



III. WKB EXPANSION FOR SLOWLY VARYING POTENTIALS 



Let U — Uh{x) be a real valued potential on some interval, with dxU — 0{\h\) and &lU = 0(h'^), where \h\ ^ 1 is a 
parameter. Let £ be a real energy and assume that lAh{x) — £ is bounded from above and from below by some strictly 
positive constants that are independent of h. This means that we are in the classically forbidden region. Then 



[-dl +U-£\{U- £)-^'^ e- r ("-^)' 
16 

-fd.' (U-Sf- ^^^2-^ 



-^^{U- £)-"' {djuf + \{U- £)-'/' dlU 



' P dx' (U-S f 



and therefore 



(3.1) 



(3.2) 



is a good approximation to a corresponding exact eigenfunction, even over intervals of length Od/ij"^). 

In the following sections, we will apply the above approximation in the barrier regions [a, b] and [c, d] with h equal 
to the X derivative oiV + W \n these intervals. 
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IV. THE DRIVING TERM 



Equation (2.12) can be solved by evaluating the instantaneous eigenstates of the potential Vmi + W. We put 
fl{x, t, E) = eKp{—iEt) fj,(x , t, E) and suppose that AV and s are slowly varying functions of time so that also t, E) 
is slowly varying in time. Thus in the equation 



[-idt -dl + Vmix) + Wis, x) - E] fi{x, t, E) = 0, 



(4.1) 



we make a very small error if we neglect the term —idtfi, as we shall do in the following. In the emitter region x < a, 
we take fj,(x, t, E) as the sum of a left and a right-going plane wave at energy E 



fi{x,t,E) 



ly/E{x — a) 



r{E)e 



— iy/E{x — a] 



(4.2) 



where r{E) is a reflection amplitude to be computed. Note that the normalization factor in (12) is chosen in order 
to have / dx fi{x, t, E)^{x, t, E') — 5{E — E') in agreement with t he e xpression of the electron density (2^) in terms 
of an integral over the energy E. We propagate the expression (4.2) to the adjacent regions by requiring /i to be 
of class and applying the WKB approximation described in Section III. In the interval [a, b] the potential is 
Vau + W = Vo + a{x — a), where 



_ Snag^s (d-c)~ AV 
b — a + d — c 

plays the role of the small parameter h of Section III. For a < a; < 6 we can then use the WKB approximation 

1 {Vo~Ey/^ ^^^^ (Vo+c.ix'-a)~Ey^' 



fj.{x,t,E) = 



t{E) e 



^i^Tl (Fo + a{x ~a)- i?)i/4 
where t{E) is a transmission amplitude to be determined with r{E) from the condition at a; = a 



1 + r{E) = t{E) 
iVs - iVs r{E) = t{E) 



{V^~Ef'^--^{V^-E)-^/\ 



Neglecting the last term in the bracket, which is Odaj), we get 



r{E) 
t{E) 



1 -^(vb/i;- 1)1/2 

2 

l-i{Vo/E-lY/^' 



(4.3) 
(4.4) 

(4.5a) 
(4.5b) 

(4.6a) 
(4.6b) 



Note that the neglected term would give correction factors 1 + C(|a|) to r{E) and t{E). 

At X = b we can set up a similar transition problem but here Vmi + is continuous and the corresponding 
transmission amplitude is 1 + Odaj). Neglecting again a factor l + Odal), for6<x<cwe get 



where 



^ioit,E) = 



fi{x,t,E)=fio{t,E) e-(Vo+a(6~a)-£)i/^(.^b) 

(Vb - £:)l/4 2 Q[iVo-Ef/^~{Vo+a{b-a)-Ef/^]2/3a 



yi;^ (Fo + a(a; - a) - S)i/4 1 + i{Vo/ E - 1)^/^ 

Only this expression of fj, in the region [6, c] will be used in the following as driving term of Eq. ( ^.13 ). 



(4.7) 



(4.8) 
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V. RESONANCE AND RESONANT STATE 



In this section we will obtain a WKB approximate expression for the ground state resonance A(s) = Efi{s) — iT{s)/2 
and the corresponding resonant state e{s,x) for the potential V + W. We will assume that c — 6 is bounded from 
below and from above by positive constants, while b — a and d — c are large enough. 

To start with, we recall the construction of the ground state eigenvalue Eq of the potential Vw{x) — 
Vq [1]-oo,6] (2;) + l[c,+oo[(2;)] which coincides, up to the constant shift 



AE = 



Sirag^s (b ~ a){d ~ c) ~ AV{b - a) 



b — a + d — c 

with the potential V + in the well region [b, c] . The corresponding ground eigenstate is 



(5.1) 




6)/2) e-(^°-^o)'^'(b-x) 
Efix-{b + c)/2)) 

b)/2) e-^^^-^o f'i^-c) 



X <b 
b < X < c 
X > c 



(5.2) 



where < Eq < min(Vb, tt^/ (c — fe)^) is determined by the requirement that e^(a;) is of class 

tan (/E™(c - 6)/2) = [Vo/E^ - l)l/^ 



and the normalization constant is 



— 



EH 



c-b ^ {Vo-E^)^/^ 



Vo{Vo~E^y/^ 



Vo 



-1/2 



(5.3) 



(5.4) 



where we used the identities cos^ u — {1 + tan^ u) ^, (sin2u)/2 = sin tt cos it = tanu (1 + tan^ u) ^. In the following 



we will assume that EJl 



AE <Vn- AV. 



Next we look at the ground state of the potential 



Vo + a{a — 6) x < a 

Vo + a{x — 6) a < X < b 

Vb{x) ^ { b < X < c 

Vo + j3{x - c) c < X < d 

Vo + l3{d-c) X > d 



(5.5) 



whic h co incides, up to the constant shift AE^ with V + W on the larger region [a,d] which includes the barriers. In 
Eq. (5_^) a is given by (|4.3| ) and 



inag^s (b-a) - AV 
b — a + d — c 



(5.6) 



Note that the potential Vb has been obtained by bending the barriers of V^ in the intervals [a, b] and [c, d] proportionally 
to a and /3, respectively. Let Eq be the ground state of Vb and eg (a;) the corresponding eigenfunction. Since 
\a\ and \(3\ are small, from the same WKB considerations of Section III we have Eq — Eq + 0{\a\ + |/3|) and 
6o(^) = ^oi^) + ^(l*^! + To get the leading asymptotics of the resonance width, we need to determine the linear 
contribution to 0{\a\ + \(3\) in Eq. By differentiating the eigenvalue equation for the potential Vb, we have 



-l-oo 



dx eg (a:) daVb{x)\^^^^QeQ{x) 



dx ix-b)\e^ix)\' 



(5.7) 



9pE'q\^=P=q= dx elix) dpVb{x)\^^^^Qe'Q{x) ^ dx {x - c) \e}^ {x)\' (5.8) 



and using ( |5.2| ) we get 
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a=/3=0 



J —OO 



Wo{Vo - 



Observing that a — (3 — STra^^s, we finally get 



(5.9) 



(5.10) 



The real part i?/?(s) of the shape resonance of —d^ + V + W which is close to the ground state eigenvalue of 
—d'^ + Vb + AE is very well approximated by the above calculated Eq + AE which can be rewritten as 



where 



and 



V = 



Er{s) = Er{0) + r,s, 
ER{Q)^E^-AV{b^a)/e 
87:ag\b - a){d ~ c) %T:ag'{C^' fE^ 



b — a + d — c 



WoiVo - E^) 



(5.11) 
(5.12) 

(5.13) 



Now we discuss the determination of the imaginary part r(s) of the resonance. In the interval [a, d] the ground 
state of Vb is 



' cos(/Eg(c-b)/2)(Vo-j;g)V^ - f dx' {Vo+a{x'-b)-El)^'^ 
(Vo+a(x-h)-E^^Y/'i ^ " 

el{x) = { COS (/E^(a; - {b + c)/2)) 

cos(/Eg(c-b)/2)(yo-i;g)i/- ^ dx' {Vo+l}{x'-c)-E''^)^'^ 
iVo+Pix~c)-E^y/i ^ " 



a < X < b 
b < X < c 

c < X < d 



(5.14) 



where Cg — Cq + 0{\a\ + \(3\). In the interval [a,d], the resonant state e{s,x) can be approximated by adding to 
(5.14) terms due to reflections at x = a and x = d. For x ^ d we try with 



e(s, x) 



cos (/E^(c - b)/2) {Vo - Elfl^ _ J. 



{V„+fJ{x'-c)-E^„)'/^ 



{Vo + p{x - c) - E'^y/^ 



X e 



l^^-~{Vo+l3{d^c)~E'g)^/^{x-d) _^ _^^{Vo+l3{d-c)-E^„y^^{x~d) 



(5.15) 



where we have also replaced the exponent with its linear approximation aX x — d. For x d we try the right-going 
plane wave 



Cg'cos(7^(c-6)/2) {Vo-E^of/^ 



■c)-^^o)''% ^^(E1,-P{d-c)f'Hx-d)^ 



The condition at x = d gives, up to terms 0(1/31), 

l + r = t 

^{Vo + (3{d - C) - 4)1/2 ^ ^y^ ^ _ ^) _ ^b)l/2 ^ ^ .(^6 _ _ ^))l/2 ^ 

which determines r and t so that for x }^ d we have 



(5.16) 



(5.17a) 
(5.17b) 



e(s, x) 



X exp 



Co- cos (/^(c - 6)/2) (t/o - i?o')'/' 



(K) + /3(d-c)-i?o^)i/4 



. (i^o^-/?(d-c))V2 
{Vo+f3id-c)-E',y/^ 



2 r 

3^5 



(K, - i?o')'/' - (^0 + P{d - C) - 4)3/2] ^ ^(^6 _ ^(^ _ ^))l/2(^ _ ^) 



(5.18) 
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In these calculations we have assumed that — P{d — c) > 0, Vq + (3{d — c) — Eq > 0. The first inequality is always 
fulfilled in experimentally relevant situat ions, while the second one, equivalent to Efi(s) < Vq — AV^ may be more 
critical and, possibly, one should replace ( 5.181 ) by a more complicated formula. 
The same calculation can be repeated for x = a. For x ^ a we try with 



X |^g(Vo+a(a-6)-£;S)i/2(^_„) ^ ^^-(Vo+a{a-b)-E'^)^/^{x-a)'j 

with a new reflection amplitude r. For a; < a we try the left-going plane wave 



(5.19) 



e(s, x) 



C^cos 



(v/^(c-b)/2) {Vo-E'„y/^ ^_^.^^, 



(Fo + a(a-6)-i?o^)i/4 
with a new transmission amplitude t. The condition at a; = a gives, up to terms Odaj), 

1 + r = t 

{Vo + a{a - 6) - E'^y/^ - (Vb + a(a - b) - Elfl"" r = -i{E^Q - a{a - b))^'^ t 
which determines r and t so that for x ^ a we have 



(5.20) 



(5.21a) 
(5.21b) 



cos (/^(c - b)/2) [V^ - ElYl^ 
e{s,x) = , ^ ^u..,. 2 



(Fo-a(6- a) -^0^)1/4 



1 - i 



. {El + a{b-a)f/^ 



{V^-a{b-a)-ElY/^ 



X exp 



2 r 

3a 



(Fo - a{b - a) - E',f'^ - (Fo - E'^f'^] - i{eI + a{b - a)Y/-\x - a) 



(5.22) 



Note that for x < a, e(s ,x) i s a true left-going plane wave only for not too large when Eq + a{b — a) > 0. If 
Eq + a{b — a) < 0, Eq. ( 5.22| ) becomes an exponentially decaying function whose correspon ding p robability current 



ing p r 

density vanishes. Since Eq + a{b — a) = Efi{s), this case corresponds to Efi{s) < 0. In Eq. ( 5.22| ), we also assumed 
that Vo - a{b - a) - E^ > 0, i.e., Eh{s) < Vq. 

The resonance width can be now computed by means of the Green formula 



r(s) / dx\e{s,x)\'^ — 2lra(^e{s,x)dxe{s,x) 



d' 



(5.23) 



where a' < a and d' > d. The integral in the l.h.s. of ( 5.23| ) is 1 + 0(101 -I- |/3|) and using ( 5.1^ ) and ( 5.22| ) we get, up 
to such a factor, 

T{s) = %{C^f eUVo-EIY'^Vq-^ 

(Vq + p^d - c) - ElY'^ {El - p{d - c))i/2 ^[iVo-Elf''-[Vo+m-c)-Elf'^]A,^P 



+ {Vo ~ a{b ~a)~ Elf/^ {El + a{b ~ a))f e[iVo-^i^-a)-E^,?'^-(Vo~Elf'^Y/Zc. 



(5.24) 



where we used u+ = 0{u) u. 



VI. ONE MODE APPROXIMATION 



Equation (2.13) can be simplified by developing v into the instantaneous eigenstates of the potential V + W 
and keeping only the contributions from the discrete resonant states, i.e., neglecting the contributions from the 
continuous spectrum p7| . For the moment, we will suppose there is only one resonant state and put v{x,t^E) = 
exp{—iEt)z{t, E)e{s, xjwhere e(s, x) is the (ground) resonant state of the potential V + W 
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[-A(s) ~dl^ V{x) + W{s, x)] e(s, x) = Q 



(6.1) 



with complex eigenvalue A(s) = Eii{s) — iT{s)/2. The eigenfunction e(s, x) is of class on the contour 7 = 
(e*^] — 00, 0] + a) y [a, c?] U (d + e'^[0, +00 [) for 6 conveniently chosen and satisfies 

(6.2) 



dx e{s,x) =1, 1 dx e{s,x) dse{s,x) — 0. 

7 -^7 



Multiplying ( 2.13 ) with e{s,x) and integrating over 7, we get 

dtz{t, E)^i[E~ A(s)] z{t, E) + S(t, s, E) 

with the driving term given by 



B{t,s,E)~iVQ I dx ii{x,t,E)e{s,x) 

Jb 



and the sheet density (2/7) reduced, with small error, to 

s(t)= / dEg{E) \z{t,E)\^ = \\z{t)f. 



(6.3) 



(6.4) 



(6.5) 



VII. FIXED POINTS AND LINEARIZATIONS: GENERAL RESULTS 

We consider the vector field in the l.h.s. of (E^), 



where 



Viz,E)^Ai\\zf,E)z{E)+B{\\zf,E), 



A{s, E) - -r(5)/2 + i{E- {Er{Q) + r,s)) , 



(7.1) 



(7.2) 



is a non-vanishing function. For simplicity, we assume that B is independent of t. When B varies slowly with t, the 
discussion below should be applied to each such fixed value of t. 

We first look for fixed points of V, i.e., functions z = z{E) in L'^{g{E)dE) with V{z{E),E) = 0. Clearly z = z{E) 
is a fixed point iff 



z{E) 



B{\\zr,E) 

A{\\z\\\Ey 



so the L^-norm s ~ has to satisfy 



s=,dEgiE) Mf:^. 



Conversely, if s > is a solution of ( |7.4D , then 



z{E) = 



B{s,E) 
A{s,E) 



(7.3) 



(7.4) 



(7.5) 



gives the unique fixed point of V with ||z|p — s. 

Assuming that we have found a fixed point z = z{E), we look for the linearization of the vector field V at that 
point. By giving an infinitesimal increment 6z{E) to z{E), the corresponding increment SV to V is 



SViz, E) = A{s, E)5z[E) + {{5z\z) + {5z\z)) {dsA{s, E)z{E) + dsB{s, E)) , 



where s — is the corresponding solution of (7.4) and {u\v) = j dE g{E) u{E)v{E). Hence, 
5V(z, E) = A{s, E)Jl{E) + {{5z\z) + (fc|z)) {dlA{s, E) z{E) + 'dlB(s, E)) , 



(7.6) 



(7.7) 
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so with u{E) = 5z{E) and v{E) — 6z{E), we get the complexification of the hnearization, 



A 0_ 
A 



{{u\z) + {v\z)){chAz + dsB) 
{{u\z) + {v\z)){d,Az + dsB) 



(7.8) 



The matrix in the first term of the r.h.s. has continuous spectrum contained in —T{s)/2 + iR and the second term 
appears as a rank one perturbation. If A e C is an eigenvalue of C with real part different from —T{s)/2, we get 



{A{s, E) - \)u + {{u\z) + {v\z)){d,A{s, E) z + dsB{s, E)) = 



(7.9a) 



{A{s, E) - \)v + {{u\z) + {v\z)){dsA{s, E)z + dsB{s, E)) = 0. 



We must then have 



u{E) = K 



dsA{s,E) z + dsB{s,E) 
A{s,E)-\ 



(7.9b) 



(7.10a) 



i{E) K 



dsA{s,E) z + dsB[s,E) 



A{s,E)~X 

{v\'z). In order to have a non-trivial solution k 7^ 0, it is necessary and sufficient that 

^z. ds{{A-\)(A-\))\B\^ -A{A-\)BdsB-A{A-X)BdsB ^ 
cLE g{E) J r — 0. 

{ReA- \f + {lmAf) \Af 



where k — {u\z) 
1-f 



Here, the l.h.s. is real for real A, and tends to 1, when A +00. 

On the other hand, the s-derivative of the l.h.s. minus the r.h.s. of IS 



/ 



dE g{E) 



\B{s,E)\^d, \A{s,E)\' - \A{s,E)\'ds \B{s,E)\^ 
\A{s,Et 



(7.10b) 



(7.11) 



(7.12) 



which coincides with the l.h.s. of (7.11) for A = 0. So if the expression (7.12) is < 0, we see that (7.11) must have 
a solution A > 0. Let us say that the fixed point is stable if the spectrum of the linearization C is contained in the 
open left half-plane and unstable otherwise. The discussion above then gives: 



Proposition VI I. 1 Let z be a fixed point of V so that (7.4) and ( [7.4 /) hold. If the s-derivative of the l.h.s. minus 
the r.h.s. of (7^) is < 0, then z is an unstable fixed point. More precisely, the linearization C then has an eigenvalue 
which is > 0. 



VIII. FIXED POINTS AND LINEARIZATIONS: THE SMALL-F LIMIT 



In this section we assume that the driving term B{s, E) is a sufficiently smooth function of E, at least near the 
point Efi(0) + r]s, where s solves (7.4). When the barriers are very wide, r(s) will be very small and 



1 



1 



\Ais, E)\^ im/2f + {En{0) +r^s- Ef 



is a function of E which is sharply peaked at Ejj{s) = Efj{0) + ijs. In (|7.4|) it is therefore justified to replace 
g{E)\B{s, E)\'^ by the constant value g{ER{0) + ■ns)\B{s, Er{0) + r]s)\'^. Then (|7.4D is well approximated by 



s = 2t: 



gjEnis)) \B{s,En{s))\' 
Vis) 



(8.1) 



We shall next apply a similar argument to the equation ( [7.11 ) for the eigenvalues of the linearization C and, for 
more transparency, we start with a simplified case, in which 
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B and F are independent of s. (8-2) 



In this case, (7.11) reduces to 



, . / ,p {E~En{0)-Vs)g{E)\B{E)\^ ^ 

^1 [{T/2 + Xr + {E-En{0)-7^smT/2r + {E-ER{0)-rjsr] ^""-^^ 



We shall use, 

r'+oo 



t2 r Rep > 0,Reg > , , 

(<Z2 +t2)(p2 +^2) I Rep > 0,Reg < ^^'^^ 

+°° 1 f ,\ Rep > 0,Reg > 

(g2 +<2)(p2 +^2) I Rep > 0,Reg < 0. ^^'^^ 



If we replace g{E)\B{E)\'^ in the integral in ( p.3[ ) with its value at i? = £"^(0) + 77s, that integral vanishes si nce the 



integrand becomes an odd function of £^ — Efi{0) — rjs. Instead, we get an approximation of t he integral in (8.3) by 
replacing g{E)\B{E)\'^ with the hnear term in its Taylor expansion at i? = i?fl:(0) + 77S. Using (^.4|), we then get from 

(ED 

, _ .^milEmiJE^ = 0, when r/2 + ReA > 0, (8.6) 

1 + A 



, ^ 2, (^W(Mto)l . 0, when r/2 + ReA < 0, (8.7) 
A 



where {g\BY)' = dE{g\BY). The solution of (|J|) is 

A = 27rr,(5|6n'(i?^(0) + 7,5) - T, (8.8) 
and this is an eigenvalue of the linearization C as long as 

2nT^{g\B\^)'{Enm+Vs)>^. (8.9) 



The solution of ( |8.7| ) is 

A = ~2^77(5|6n'(i?fl(0)+77s), (8.10) 



and describes an eigenvalue of C precisely when ( p.9[ ) is fulfilled. We then have the following conclusion under the 
simplifying assumption (|]^) and in the small-F limit. 

When 27r?7(g|Sp)'(i?fl(0) + 77s) < r/2: no eigenvalues of C and hence an attractive fixed point. 
When r/2 < 2TTri{g\B\'^)' {Eii{0) + 7/s) < F: two eigenvalues of C and still an attractive fixed point. 
When 2'Krj{g\B\'^)' {E + 77s) > F: two eigenvalues and a non-attractive fixed point. 
The main conclusion under the same assumptions is then: 

Proposition VI II. 1 We have an attractive fixed point precisely when the s-derivative of the difference of the l.h.s. 
and the r.h.s. in (8.I) is > 0. 



Now we drop the simplifying assumption ( |8.2| ) and see that the preceding proposition still holds in the small-F limit. 
Let z be a fixed point, so that s = ||z|p (approximately) solves (B.l). In view of ([z^), Eq. (7.11) can be written as 

. , /• ,p d.aA^X)(A-X))\B\^-A(A-X)d,BB-A{A-X)Bd,^ 
^ J > [{T{s)/2 + A)2 + {E- Er{Q) - 77s)2] [(F(s)/2)2 + {E - Er{0) - 775)^] ^ ' ' 

Here, the numerator of the integrand can be simplified to 

[9,F(F/2 + X)\B\' - (F/2)(F/2 + X)ds\B\^] 

+ [{E - ErXO) - 77s)(-277|6|2 + tXid,BB-BdM - [{E - Er{0) - 7^sfd,\B\% (8.12) 
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Accordingly, we split the integral into three pieces and apply the small-F approximation to each one. If we assume, 
for simplicity, that r/2 + ReA > (which is necessarily the case if the eigenvalue A is to ruin attractiveness) we get 



1 



[dsT{T/2 + A)|gp - (r/2)(r/2 + X)dM^9 

(r/2)(r/2 + A)(r + A)_ 







r + A 

a± E ~ Eji{0) + rjs. This can be rewritten as 

27r(6>,r/r)g|gp ^TTdsjgm ^TTrjdEjgm tX7TdE{g{dsBB - BdM 

r+A r+A r+A r+A 

again at i? = £'_r(0) + rjs. Noticing that 

En{0) + vs)) = [i^dE + d,){g\B\^){s, En{0) + ^s), 
and multiplying with r + A, we get the following approximation of ( [7.11] ) 



(8.13) 



(8.14) 



X[l + tTTdE{g{dsBB - BdsB))] 

{gmis,EniO) + vs) 



-r(s) - 2tt- 



r{s) 



dsTis) + 27:-i{g\B\^){s, ERiO) + rjs)). 



(8.15) 



We assu me th at 1 + indEigidgBB — BdsB)) > 0, so that the solution A of ( S.15 ) is real and has the same sign as the 
r.h.s. of 



.15) 



On the other hand, the s-derivative of the l.h.s. minus the r.h.s. of (8.1) is 



1 - ^^{{g\B\^){s,EniO) + rys)) + 27rig\B\%s, E^iO) + r^s)^'^^'^ 



r(s) ds 

1 



r(.)2 



+ 2.^((,|Sn(., ERiO) + ,.)) - 2.M!Kfl|p±^0.r(.) 



which is of the opposite sign to the r.h.s. in (8.15). We then have: 



Proposition VIII. 2 Under the weaker assumptions above and in the small-T limit, we still have an attractive fixed 
point precisely when the s-derivative of the l.h.s. minus the r.h.s. of (8.1) is > 0. 



IX. A SIMPLIFIED DIFFERENTIAL EQUATION FOR THE SHEET DENSITY 

Consider the differential equation ( |6.3| ) 

dtz{t, E) = [-r(s(0)/2 + i{E - (ErXO) + 7j.sit)))] z{t, E) + S(s(t), S), (9.1) 

where s{t) ~ \\z{t, •)|p, and where we could also let S be a slowly varying function of t through s{t). Assuming s{t) 
to be a known function, the solution of (9.1) with a prescribed initial value at time t^ is 



J to 

t{E~EniO))(t~to)~ r dt'r(s(t'))/2~iv r dt's(t'} 

+e -^'0 z(to,Ej. 

Assuming that the solution has existed as a bounded solution for a very long time, say from the time — c», we can let 
to tend to — oo in the formula above and get 



z{t,E) 



dt'e 



t(E-ER{0)){t-t')-f^^ dt"T{s{t"))/2-ir, J*^ E) 



(9.2) 
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Taking the scalar product of (9.1) and z gives the following equation for the derivative of the sheet density 



where 



—s(t) = 2Rc(z|9tz) = -r(s(t))s(i) + 2Re(z|S), 



2Re(z|6) = 



(9.3) 



(9.4) 



We now assume that s{t) varies slowly with t and replace B{s(t'), E) in the above integral by S(s(t), E). Making the 
^^-integration first, we get 

2Re{z\B) = 2Re / * dt'e"'"'"^")^*-*''-/*' ^*"r(-(*"))/2-"; J/, f _ t)^ (9.5) 

where denotes the Fourier transform with respect to E. Assuming g{E)\B{s{t),E)\'^ sufficiently smooth as a 
function of E, we see that !F{g\B\'^){s{t), t' — t) is sufficiently rapidly decreasing as a function of t' — t for the following 

approximations to be made: i) since T{s) is small, we may assume that exp|— J^* dt"T(t")/2^ ~ 1 ii) since s{t") 
varies slowly, we may replace J^* dt"s{t") by s{t){t — t'). Wc then get 



2Re(z|S) ~ 2Re / dt'e-^(^«(°)+''"(*»(*-*').F(5|6n(s(t), t' - t) 

J — OO 

= 2Re / dt'e*(^«(")+''"(*»*^(ff|Sn(s(i),t')- 



Using the property !F{u){—t) = T{u){t), valid for any real valued function u{E), we obtain 



+ 00 



2Re(s|6) = / dt'e^(^«(°H'7^(t))t>(g|s|2)(,(i)^ = 2^{g\B\'){s{t), Er{0) + ^s{t)). 



Inserting this in ( |9.3D , we get the approximate differential equation for the sheet density s{t) — \\z{t, •)|| 



j,sit) = -r(.(0) 



s{t) - 2ti 



{g\Br){s{t),Ej,{Q) + ^s{t)) 



Visit)) 

This equation is valid for slowly varying solutions which have evolved for a time much longer than T^^. 



(9.6) 



■7) 



X. QUALITATIVE DISCUSSION AND NUMERICAL RESULTS 

We start by examining the simplified fixed point equation ( ^.ij ). For < i? < Ep with Ep <C Vb, we have 
Vq — E ^ Vq (of the same order of magnitude). By evaluating the integral in ( ^ ) with e(s, x) approximated by ( |5.2| ) 
and the driving term given by (4.7), we have 



Assuming for simplicity zero temperature, so that giE) — 9iE)iEp — E)+/2'k, we get 

giERis))\Bis,ERis))\^ 
~ iC^)^V-'E}^'ERis)^_l\Ep - i^^(s))+e[(^«-^-(^))'''-(^»-^«(^)+"(''-'^))''l4/3o 



Recalling that Enis) = i?j?(0) + rjs = + aib — a) ^ Eq — /3(c? — c) — AV, where a and /3 are given by (4.3) and 
(^.6D, respectively, from (5.24) we get 
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We will consider the following two cases: 

1) The barrier [c, d] is more opaque than [a, h] in the sense that the exponential factor in the second term of the above 
expression for T{s) is much smaller than the exponential factor in the first term. 

2) The barrier [a, b] is more opaque than [c, d]. 

In the intermediate case when the two barriers have opacity of the same order, the discussion of case 1) will roughly 
apply. Notice that opacity depends not only on the relative sizes of 6 — a and d — c, but also on s and AF. Therefore, 
we may have transitions between the two cases when these parameters vary. Interesting phenomena appear when the 
case 1) is possible and we start with that case, recalling that Eji{s) = £'^(0) + rjs = Eq — AT^(6 — a)/£ + 'qs. In this 
case (and neglecting, to start with, the possibility of a transition to the case 2)) the first term in the expression for 
r(s) dominates, except when Er^s) is negative or very small and positive. The function 



f{s) = 2^ 



g{ER{s))\B{s,En{s))[^ 
T{s) 



vanishes for Er{s) < and rises very sharply (with a square root singularity at Er{s) = 0) from to 

/max ~ Ep 



(10.1) 



(10.2) 



when Eji{s) is increased from to a small positive value. When Efj{s) is further increased, the function /(s) decreases 
roughly linearly and vanishes for Eii{s) > Ep. The values Eii{s) — 0, Er{s) = Ep correspond to 



s={AV{b-a)/e-E^)/^, 



[AVib - a)/i - i?- + Ep)/r^, 



(10.3) 



and describe the boundary points of the support of the function ( |10.1D . When AV is increased, these two points 
move to the right with the same speed as shown in the example of Fig. 2. In Fig. 2 we also see the graphical 
solution of Eq. (B.l), s — /(s), for different values of AV. It is clear that ( ^.l[ ) will first have only one solution when 



AV{b — a)/£ — Eq < 0, then three solutions for AV in some interval, until AV(b — a)/£ — Eq ~ ?7/max, and again 



only o ne so lution for even larger values of AV. According to the results of section [VIII , we see that in the case in 
which ( ^.l|) has only one solution, this solution corresponds to an attractive fixed point, and when there are three 
solutions, the smallest and the largest of these correspond to attractive fixed points, while the intermediate solution 
corresponds to an unstable fixed point. 

For many experimentally relevant situations the resonance width is much sma ller than the other energy scales 
(essentially Ep). In this case we m ay expect the simplified fixed point equation (3.1) to be a very good appr o ximat ion 
of the more correct equation (7^), except near the boundary points of the support of t he f unctio n (10.1 ). This is 
confirmed by Fig. 3 where we show the numerical solutions (st able and unstable) of both (7.4) and (3.1) for a system 
having r{0)/ Ep ~ 0.01 at AV — 0.2 eV. In the case of Eq. (7.4), the corresponding energy integral has been evaluated 
on a uniform energy mesh having a density of points 3> r(0)^^. 

The solutions of the simplified differential equation ( |9.7|) converge to one of the solutions of (8T), associated to 
an attractive fixed point. The phenomenon of hysteresis then becomes clear. We let AV increase very slowly from 
some sufficiently small value up to some sufficiently large positive value and subsequently decrease it very slowly, back 
to its initial value. Consider a corresponding solution of the time dependent Schrodinger equation ( ^.3|) so that we 
expect the corresponding evolution of the sheet density to be well approximated by the solution of ( |9.7| ), where B 
varies slowly with time. First, there is only one (attractive) fixed point and the time dependent solution has to stay 
close to that fixed point. Then we have creation of a pair of fixed points (one stable and one unstable) away from 
the solution, but the solution continues to stay close to the old (stable) fixed point. When AV reaches a sufficiently 
large value, the unstable fixed point runs into the old stable one and they both disappear. At this point, the time 
dependent solution has no other choice than to converge to the only remaining fixed point (which is stable). When 
decreasing AV back to its initial value, we have the same behaviour, in the sense that the solution stays close to 
the initially unique fixed point as long as it exists and converges to the new unique fixed point after the old one has 
collapsed with the unstable one. The bias energy AV at which this collapse happens depends on the value of the time 
dependent solution and therefore is different when AV is increased or decreased. 

The phenomenon of hysteresis is clearly seen in Fig. 3, where the collapse points for AV^ decreased from large 
values and AV increased from small values have been marked with A and B, respectively. We have AVa < AVb- We 
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can estimate the order of magnitude of the hysteresis width AVb - 
the condition Efi{s = 0) = and AVb by the condition Eb,{s ~ /„ 

AVb - AVa ~ vfm.J{b - a)-' ^ a 



AVa by considering that AVa is determined by 
,x) ~ 0. We have 



,'Ep{d-c). 



(10.4) 



In Figs. 4 and 5 we show the time dependent evolution of the sheet density s{t) when we start from a fixed point 
solution corresponding to the point A ot B and give an instantaneous small decrement or increment SV to A Va o r 
AVb, respectively. In these figures, the thick Hnes are the solutions of the full Schrodinger equation (3^) and ( |6.5| ) 
and the thin lines the solution of the simplified differential equation (5/7). In Fig. 4 the solutions corresponding to 
the small-F limit and the full Schrodinger equation start, as shown in Fig. 3, from different fixed point values, s(0), 
and converge to the same (approximatively) values. On the other hand, when the starting point is B (Fig. 5) the 
small-F approximation is close to the solution of the full equation except for the value which s{t) has to converge to, 
again in agreement with Fig. 3. 

As a third example of time evolution of the sheet densit y of electrons in the well, in Fig. 6 we show the behavior 
of s{t) solution of the full Schrodinger equation ( |6.3| ) and (6^) after an instantaneous change SV of the bias energy 
corresponding to the point C of Fig. 3 well inside the hysteresis region. If \SV\ is chosen sufficiently large, we observe 
oscillations of s{t) on the picosecond time scale. Contrary to the claim of ||lj], these oscillations are damped since 
s{t) has to converge to the fixed point solution corresponding the bias energy AVc + SV. 



In the case 2), when the barrier [a, b] is more opaque than the barrier [c, d], the function (10.1) is very small, and for 



solutions of (8J) we can observe only a microscopical hysteresis effect, du e to the square root singularity at i?_R,(s) = 0, 
which is likely to be completely absent in the more correct equation (7.4). The absence of the hysteresis effect in this 
case is in agreement with the experimental results of [|l8| and is discussed in [|l^ . 

Let us finally consider the case of very wide barriers and see that a transition between the cases 1) and 2) has to 
take place in the hysteresis region . Let c — b = constant , [b — a)/(d — c) = constant < 1 , and b— a —^ oo. In this limit , 
? 7 ^ ap^(6— g), and the values in ( 10. 3| ) are the endpoints of a short interval of length ~ EpaB / [b— a). Le t us consider 
( ^■l[ ) when AV is increased from the initial value {E}^ — EF)i/{b — a) for which the right end point in ( 10.3 ) is 0. If 
the constant {b — a)/{d — c) is sufficiently small, we are in the case 1). For (E^ — Ep)l/{b— a) < AV < EQ£/{b~ a), 



we remain in the case 1), provided that (b — a)/{d — c) is sufficiently small, and (8.1) has a unique solution. At 
AV = EQ£/{b — a) we have the creation of two new fixed points. If we follow the old fixed point, we cannot remain 
in the case 1) until it disappears. Indeed, if we did, the disappearance would take place when s ^ /max and at a bias 
energy AV ~ ag^{d — c)/max obtained by the condition Er^s ~ /max) — 0. Since Er{s) = Eq+ a{b — a) is between 
and Ep, the inclination a of the first barrier would have to be very small and we get a finite inclination /3 ^ ~o-~b /max 
for the barrier [c, d]. Therefore, when 6 — a — > cx3 only the opacity of the first barrier would diverge and, at some 
point, we would be no more in the case 1). What will actually happen is that when AV reaches some value which 
is bounded independently of 6 — a, we have a trans ition from the case 1) to the case 2) and /max decreases to some 
value which is much smaller than the r.h.s. in (10.2). This will cause the disappearance of the fixed point for a much 
smaller value of s. When a transition from case 1) to case 2) happens, we still observe a hysteresis phenomenon, but 
this is now caused not only by the translation of /(s) as a function of AV but also by a variation of its height. This 
effect is already apparent in Fig. 3 where we see a decreasing of the height of /(s) when increasing AV . 



XI. THE CASE OF SEVERAL RESONANCES 



In this section we discuss very briefly the case with several shape resonances. Much of the discussion is similar to 
the case of one resonance and we shall assume that we are in a parameter range where all the WKB considerations 
apply. 

First we review the approximation for the shape resonances. We start with the potential Vw and consider its 
eigenstates e™(x), j = 0, 1, . . . , — 1 and the corresponding eigenvalues < Eq < E]" < . . . < E^_^ < Vq. Since 
ej'{x) is even as a function of x — {b + c)/2 for even j and odd for odd j, we have 

f sin((j + l)7r/2 - /Ef(c - b)/2)e-^^°-^T)''"i''-x) ^ ^ ^ 

ej{x) = Cf I sin((j + l)7r/2 + /E™(x - {b + c)/2)) b<x<c 
[ sin((j + l)7r/2 + /Ef (c - b)/2)e-^^°'^T^'"i^-''^ x > c. 

The condition at a; = 6, or equivalently at x = c, gives the quantization condition 

tan ((j + l)7r/2 + ^{c - 6)/2) = -{Vo/EJ - 1)-'/^ 
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which can also be written as 



tan I 



(^(c - b)/2 - jV2) - [Vo/EJ - 1)1/2. 



(11.1) 



Representing this equation graphically, we see that — 1 is the largest integer > 1 with \/Vq{c — b)/2 > {N — 1)tt/2. 
The functions ej{x) are normahzed, if we choose 



-1/2 



^ \Vo{Vo-Efy/^ 



Vo 



(11.2) 



The eigenvalues E^ associated to the potential Vf, in {p.5\) can be studied as before, and we get 



^ Wo{Vo-E^ 



(11.3) 



In the following, we neglect the error 0{o? + 0^). The shape resonances Aj(s) = Erj{s) — iTj{s)/2 for the potential 
V + W are then given by 



where 



and 



ER,,iO)^Ef -AVib-a)/£, rj, 



8TTa^\b - a){d ~ c) ^nag\CffEf 
b-a + d-c Wo{Vo-Ef) ' 



(11.4) 



(11.5) 



T,{s) = 8{Cfr E^ {Vo-E^y/^V,-' 

{Vo + Pid - C) - E'^y/^E'^ ~ Pid - c))l/2g[(yo-iJ,^)^/^-(yo+/5(d-c)-i5^^)3/^]4/3/3 

~{Vo ~ a{b -a)- E]Y'\E] + a{b - a))^/2^[(^''-"(''-)-4)-^''-(^o-B,^)^/^]4/3o 



(11.6) 



The corresponding resonant state ej(s,a;), satisfying (i^), can be described as in section 0. 



We still try to solve (2.11) in two steps. Equation (2.12) is treated as before, while the Eq. (2.13) is now handled 
by letting V he a, linear combination of the N resonant states eo(s, x), . . . , eAr_i(s, x). More precisely, we write 
i'{x,t,E) = exp{—iEt)i>{x,t,E) and 'jl{x,t,E) = exp{—iEt)ii{x,t, E), so that ( p. 13 ) becomes 



^idt -dl + V{x) + W{s, x) - E]i^{x, t, E) ^ Va l[b^c] t, E). 



Assume, 



'{x,t,E) ^ ^ Zk{t,E)ek{s,x), 



(11.7) 



(11.8) 



A;=0 



where s is defined in (2.7) and hence will be time dependent. The functions eo(s, x), . . . , ejv-i(s, x) approximately form 
an orthonormal family in L'^{[{a + b)/2, {c + d)/2]), and if we assume that v dominates over /i in [(a + &)/2, (c + d)/2] 
then, with a small error, we have 



7V-1 



,(o = Eii^fc(*'-)ii' = ii^(^'-)f: 



(11.9) 



fc=0 



where the norms are in L'^(g{E)dE) and in L'^{g{E)dE)^ , respectively. 

Substituting (11.8) into ( 11.7] ), multiplying by ej{s,x) and integrating over the contour 7, we get 
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N-1 ^ 

dx [-idt + Xk{s) - E]{zk{t,E) ek{s,x)) ej{s,x) ^Vo / dx fj,{x,t, E)ej{s, x). 
k=o ''^ •'^ 



(11.10) 



From the relations dx efe(s, x)ej (s,x) — 6k. j , we conclude that dx {dgBkis, x))ej{s, x) is an anti-symmetric matrix, 
and since ek{s,x) are approximately real functions near [&, c], this matrix is also very close to a real one. Equation 
(11.10) can be written as 



N-l „ 

■idt + Aj(s) - E]zj{t,E) - idtis{t)) ^ / {dsek{s,x))ej{s,x) 



Va I dxii{x,t, E)ej{s,x). 



(11.11) 



Due to the facts that i) dts{t) can be expected to be very small and ii) efe(s, x) is roughly independent of s near [b, c] 
so that the integral dx{dsek{s,x))ej{s,x) can be expected to be very small, we will neglect the sum in the l.h.s. of 



( 11.11 ). In this case, we have 

dtZj{t,E) - [-Tj{s)/2 + i{E - ER^j{s))]zj{t,E) +B,{t,s,E), 



(11.12) 



where Bj{t, s, E) — iVa dx ii{x, t, E)ej{s, x). 

We assume that Bj va ry slo wly with t, so it is meaningful to look at instantaneous fixed points of the vector field 
defined by the r.h.s. of (|ll.l2 ) in L^{g{E)dE)^ . Assuming, for simplicity, that Bj are independent of t we see that 
z{E) — {zq{E), . . . , zn-i{E)) is a fixed point precisely when: 



-B,{s,E) 



-rj{s)/2 + i{E-EnAs))' 
from which we get the compatibility condition for s — WzW^ 



j = 0,...,7V-l, 



E 



dE 



9{E)\B,{s,E)\^^ 



= 0. 



(11.13) 



(11.14) 



Conversely, if s is a solution of (1L14 ), th en (11.13) defines the unique fixed point with 
In the small-F limit, as in section |VIII we get the simplified fixed point equation 



^\ j9\B,?){s,EnAs)) 



3=0 



r,(s) 



In view of (11.4), the term of index j in (11.15) is a function of s with support in the interval 



{AV{b - a)/£ - E^)/rjj < s < {AV{b - a)/£ - Ef + Ef)/!]^, 



(11.15) 



(11.16) 



and when AV^ increases this interval moves to the right with speed (6 — a)/{£rij) as sh own in the example of Fig. 7. 



In Fig. 8 we compare the corresponding fixed point solutions obtained by solving (11.14) with those obtained in 
the small-F limit (11.15) as a function of the bias energy AV. Between the points marked as A and B we observe five 



fixed points. Below we give some results about the nature of fixed points, which are more complicated than in the 
case of a single resonance and it is not clear that those results are applicable in the situation of Fig. 8. If we assume 
that they are applicable, then three fixed points are stable and two unstable. The existence of more t han th ree fixed 
points, i.e., the maximum number allowed for = 1, is related to the possibility that the intervals ( 11.16 ) are not 
disjoint, as clearly understood by Fig. 7. 

It is interesting to study the evolution of the sheet density s(t) away from a point like B in Fig. 8 where a 
(presumably) stable fixed point and an unstable one collapse while two other fixed points survive. In Fig. 9 we show 
the behavior of s{t) obtained by numerically integrating Eq. ( 11.12| ) after an instantaneous increase 6V of the initial 
bias AVb- If the total bias AVb + SV < AVc, where C is the next point where a new couple of fixed points collapse, 
s{t) converges to the fixed point closest to its initial value s(0). When AVb + SV > AVc, s{t) first approaches the 
value corresponding to the collapse point C but finally has to converge to the lower unique fixed point corresponding 
to the chosen bias. 
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Next we study the linearization of the vector field defined by the r.h.s. of ( 11.121 ) at a fixed point under the following 
simplifying assumption: 



Then (11.12) becomes 



Tj is independent of s, Bj = Bj{E) is independent of t, s, 
r]j = r] is independent of j. 



dtz,{t, E) = [-r,/2 + i{E - £;hj(0) - vs)]z,{t, E) + Bj{E). 



(11.17) 



(11.18) 



The same calculations as in s ection VIII show that the complexification C of the linearization of the vector field 



defined by the r.h.s. of ( ll.lq ) at a fixed point, is given by 

/ lio \ / [-To/2 + i{E-ER,a{Q)-r^s)]uo~ir]{{u\z) + {v\z))za \ 



C 



UN -I 



[-rAr_i/2 + i{E - Er^n-i{^) - i]s)\uN-i - iri{{u\z) + {v\z))zn-i 
[-ro/2 - i{E - En^oiQ) - t^s)]vo + n^{{u\z) + {v\z))zo 



\ VN-i ) \ [-rAr_i/2 - i{E - Eji^M^ii'^) - -qsy^VM-i + iri{{u\z) + {v\z))zn^i ) 



(11.19) 



Here, {u\z) — X]j=o {"^A^j) L'^(g(E)dE)- The operator £ is a rank one perturbation of an operator with essential 
spectrum contained in U^^^(— rj/2 + iR). Wc look for eigenvalues A G C with Re A 7^ -'Tj/2 for all j. If 
(mo, . . . , itAf-i, Wo, ... , fiv-i) is a corresponding eigenvector, we get as in section VIW 



-nZj{E) 



^ -Tj/2 + i{E~ER,j{0)~r]s)~X' ' -Tj /2 - i{E ~ Erj{0) - 7]s) - X' 

where 

K = ir]{{u\z) + {v\z)). 
Using (11.13), ( 11.20| ) in (11.21), we see that A is an eigenvector precisely when 



1 - 2,7 ^ fdE 

k=0 •' 



(E - EnAO) - Vs)g{E)\Bk{E)\^ 



[(Tk/2 + A)2 + {E- ErAO) - VsrWk/2r + {E- Er^O) - V^)^] 



= 0. 



(11.20) 



(11.21) 



(11.22) 



As in the case = 1, we observe that the l.h.s. of (11.22) for A = is equal to the s-derivative of the l.h.s. of 
(11.14). Moreover, when A — *■ +00, the l.h.s. of (11.22) converges to 1, so if it is < for A = 0, it has to vanish for 
some A > 0. Hence, as in the case iV = 1, we get: 



Pr oposit ion XI. 1 Let z be a fixed point of ( 11. IS ), so that s — WzW^ solves ( 11. 14 )■ If the s-derivative of the l.h.s. 
of ( ll.lJD is < 0, then z is not an attractive fixed point. 



We now pass to the small-F limit, where (11.14) is replaced by (11.15) and we keep the simplifying assumption ([11.17[) 



Proposition XI. 2 (small-F limit) Assume that the intervals (ll.lt) are disjoint and let z be a fixed point of 
^IIJ^ . Then z is attractive precisely when the s-derivative of the l.h.s. of ((11.13^ is > 0. 

Proof. The s-derivative of the l.h.s. of ( 11.15| ) is 



N-l 



27:7j{g\B,\^yiERA0) + ris) 



(11.23) 



where (g|6jf )' 
linearization becomes as in section IVIII 



dE{g\^]\ \ On t he ot her hand, in the small-F limit, the equation (11.22) for the eigenvalues of the 
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1-2, 

k; rfc/2+RcA>0 
k; rfc/2+RcA<0 



TT{g\BkniEii^k{0)+V-s) 
Tfc + A 

7r(ff|gfcpy(^fl,fc(0)+7?g) 
A 



(11.24) 



We are only interested in the possible existence of solutions to this equation with Re A > 0, and for such A ( 11.24 ) 
reduces to 



N-l 



k=0 



nig\ Bk\'Y{EnAO) + V''^) 
Tfc + A 



= 



(11.25) 



If A is a solution, then by the condition that the intervals ( 11.16 ) are disjoint, only one term in the last sum, say for 
k = m, is 7^ 0, so that (11.25) becomes 



while the expression (11.23) becomes 



1-27; 



Tm + A 



TT{g\B„,\^y{ER^m{0)+Vs) 
r 

^ m. 



(11.26) 



(11.27) 



It is then easy to see that the solution of ( 11.26 ) has a negative real part precisely when the expression ( 11.27 ) is 
positive, and this concludes the proof of the last proposition. 



When the intervals (11.16) have non-empty intersections, the situation is more complicated, and the following 
example is an indication that the last proposition may be false. 

Example. There exist Fi > 0, oi, 02 € .R, such that 1 — (ai/Fi + a2/r2) > 0, while 1 — (ai/(ri + A) + 
12/ {^2 + A)) — for some positive A. Indeed, choose Fi = 1, ai = 2, r2 = (5 > very small, a2 = —25. Then 
1 - ai/Fi - a2/r2 = 1 > 0. If ^ << Aq << 1, we have 



1 



ai 



0-2 



-1. 



Ti + Aq r2 + Ao 

Hence 1 — (ai/ (Fi + A) + 02/(1^2 + A)) must vanish for some A between and Aq. 

As in section p^ , we can derive a simplified differential equation for (so(i), . . . , sj\f-i{t)), where Sj{t) = ||2j(i, OIP; 
so that s{t) = X]^o^ ■Sj(0- We drop the simplifying assumption (|ll.l7| ), but keep, for simplicity, the assumption that 
Bj are independent of t. Assume that z{t,E) = {zQ{t,E), . . . , ZN^i{t, E)) is a solution of ( 11.12| ) which has existed 
for a long time with a uniformly bounded norm. As in section we take the scalar product of ( 11.12 ) with Zj and 
get 



dt 



s,{t) = -Tj{s) + 2Re {z,\Bj 



Using 



Zj{t,E) 



dt'. 



and, under the assumption that s{f) is slowly varying, we get as in section IX 

2Re {z,\B,) ^ 27r{g\B,\'){s{t),ER,,{0) + rj,s{t)) 

and the simplified equations 



-s,(i) = -r,(s(i)) 



Sjit)-27r 



{g\B,\^){s{t),ER,,{Q) + 7j,s{t)) 



(11.28) 



(11.29) 



(11.30) 



(11.31) 
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for j = 0, 1, . . . , iV — 1, and s = X^jLo^ ^j- notice that the region defined by sj > for < j < iV — 1 is stable 
under the forward flo w asso ciated to the system (11.31). Mor eover, if {sq, . . . , sjv-i) is a fixed point of this system, 
then we get precisely ( |ll.l5| ). Conversely, if s is a solution of ( |11.15| ), then 



2tt 



{g\B,\^)is,En,,{0)+^,s) 

r,(s) 



(11.32) 



defines the corresponding unique fixed point solution with s — J2j 

We end t his sec tion by investigating the linearization of ( |ll.31 ) at a fixed point solution, under the simplifying 
assumption ( |ll.l7| ). An easy calculation shows that the linearization is given by 



/ 



M 



( 



^rowo + 2^775(5|eoH'(Si?,o(0) + ns) E 



fc=0 



(11.33) 



If A is an eigenvalue of A1 with Fj + A 7^ for every j, and '(wo, . . . ,WAr_i) a corresponding non-trivial eigenvector, 
we have 



r, + A 



AT-l 



Then n ecessarily the sum is 7^ 0, and by summing these iV relations, we see that A is an eigenvalue precisely when 
(111. 251) holds. We finally get: 

Propo sition XI. 3 Under the simplifying assumption ^ILlI^ ) and in the small-T limit, let z he a fixed point of 
( 11.1^) and let {sq, . . . , Sat-i) be the corresponding fixed point solution of (11.31). Then the linearizations of (11.16) 
and ( 11.31 ) at the corresponding fixed points have the sa me eig envalues in the right half plane ReX > ( given by 
( 11.2q) )~In particular, z is an attractive fixed point for (11.18) precisely when (sq, • ■ ■ , sw-i) is an attractive fixed 
point for \ll.3\ ). 
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FIG. 1. Potential V{x) representing the band profile modified by the external bias energy AV (solid line) and total potential 
V{x) + W{s,x) including the electrostatic contribution due to electrons trapped in the well with sheet density s (dashed line). 
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FIG. 2. Graphical solution of the equation s — f{s) for different values of the bias energy AV. Note that the support of 
f{s) has width As ~ Ep/ri (equality strictly holds at zero temperature). The example shown here corresponds to a typical 
GaAs-AlGaAs heterostructure in which the parameters described in the text have the following values: nn = 2 x 10^^ cm~^, 
T = 1 K, 6 - a = 40 A, c - & = 56 A, d - c = 70 A, Vb = 0.34 eV, e = 11.44, and m* = 0.067 m, where m is the free electron 
mass. 




AV (eV) 



FIG. 3. Fixed point solutions of the sheet density of electrons in the well s as a function of the bia s energy AV^ in the same 
case of Fig. 2. The thick line is the exact case ( |7.4[ ) and the thin line the small-F approximation (8.1). Unstable solutions are 
represented by dashed lines (both thick and thin). 
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FIG. 4. Sheet density of electrons in the well s{t) as a function of time after an instantaneous decrease 6V of the bias energy 

from the point A of Fig. 3 (thick lines). The crosses are the fixed point solutions at bias AVa — 5V where s{t) is expected to 
converge. The thin lines are the corresponding small-F approximation starting from s(0) = 0. 



10 



10 



12 



11 



§ lolo 



10^ 



10"^ 



c^-^ «^ 

- 4 \ \ 

^ \ \ 

- \ \\ 

- % \\ 
~ % \\ 

\ \\ 

Is \\ 

% \\ 
% \\ 
% \ \ 






* \ * 

= \V \ 

" * \ 

' (5V=10"^ eV 

\ 


\ (5V=10 ^ eV 




\ 

Z N 

\ 






^ N 












1 1 1 1 1 1 iv 







100 



200 300 

-12 



400 



t (10^^'' s) 



FIG. 5. Sheet density of electrons in the well s{t) as a function of time after an instantaneous increase 5V of the bias energy 
from the point B of Fig. 3 (thick lines). The crosses are the fixed point solutions at bias AVb + SV where s{t) is expected 
to converge. The thin lines are the corresponding small-F approximation. For 5V not too large a ghost fixed-point solution is 
observed with s(t) decaying linearly for t <tg and tg defined by the condition Eii{s{tg)) = 0. 
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FIG. 6. Sheet density of electrons in the well s{t) as a function of time after an instantaneous change 6V of the bias energy 
from the point C of Fig. 3. The crosses are the fixed point solutions at bias AVc + 5V where s{t) is expected to converge. For 
\5V\ not too small damped oscillations are seen at the ps scale. 
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FIG. 7. Graphical solution of the equation s = f{s) for different values of the bias energy AV in a multiple-resonance case. 
We used the same parameters of Fig. 2 except 6 — a = 20A, c — b = 360 A, d — c = 50 A. 
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FIG. 8. Fixed point solutions of the sheet density of electr ons in the well s as a function of t he bias energy AV in the 
same case of Fig. 7. The thick line is the solution of Eq. (11.14) and the thin hne the small-F limit (11.15). Possibly unstable 



solutions are represented by dashed lines (both thick and thin) 
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FIG. 9. Sheet density of electrons in the well s{t) as a function of time after an instantaneous increase 5V of the bias energy 
from the point B of Fig. 8. The crosses are the fixed point solutions at bias AVb + 5V where s{t) is expected to converge. 
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